Diffusion-limited aggregation in channel geometry 
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We performed extensive numerical simulation of diffusion-limited aggregation in two dimensional 
channel geometry. Contrary to earlier claims, the measured fractal dimension D = 1.712 ± 0.002 
and its leading correction to scaling are the same as in the radial case. The average cluster, defined 
as the average conformal map, is similar but not identical to Saffman- Taylor fingers. 
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Diffusion-limited aggregation (DLA) has attracted 
considerable attention since its introduction by Witten 
and Sander in 1981 0. In this model an aggregate or 
cluster grows by capturing diffusing particles which irre- 
versibly attach to it on first contact. This is the discrete 
model of a wide variety of physical systems in the Lapla- 
cian growth class. This class can be modeled with a field, 
satisfying the Laplace equation outside of a growing clus- 
ter, where the cluster grows in proportion to the gradient 
of the field at the boundary. A generalization of DLA, 
known as the dielectric breakdown model (DBM) Q , al- 
lows growth proportional to the field gradient to the ex- 
ponent 77. DLA is regained for rj — I. The majority of our 
current knowledge about DLA is numerical (mostly in 
two dimensions and radial geometry), although progress 
has been made in the theoretical front as well (see eg. 

Si) 

One of the controversies surrounding DLA in channel 
geometry has been that the fractal dimension might differ 
from that in the radial case. This claim has been based 
on small size simulations 0, small size calcula- 
tions (S, t9j- this paper, based on extensive numerical 
simulations of off-lattice DLA in channel, we show that 
the fractal dimension is asymptotically the same in the 
two geometries. 

One of the most important differences between the two 
geometries is that for the channel the continuum version 
of the problem (Laplacian growth) has a stable solution 
without tip splitting instability and finger competition. 
These stationary translating solutions — called Saffman- 
Taylor fingers [l(j| — have been studied in viscous finger- 
ing experiments in Hele-Shaw cells ■ In this paper we 
compare them with average profiles of DLA clusters. 

Our first method for generating DBM in a channel is 
to use iterated conformal maps. The conformal map pin g 
method for the radial case is described in [l2| and |l3l| . 
For radial DLA, a map is created from the unit circle in 
the w ("mathematical") plane to the unit circle with a 
bump at a randomly chosen angle in the physical plane; 
the composition of such maps is a map from the unit cir- 



cle the DLA cluster. To produce DBM clusters, instead 
of choosing the angles at random, we use a monte-carlo 
method to select bump sites with the correct distribution 

To adapt this method to a channel, we modify the 
map of Stepanov and Levitov [l3 | for both periodic and 
reflective boundary conditions by requiring the map to 
be symmetric about the real axis. The map is given by 



(1) 



where g{w) — ^^^j^ is a map from the unit circle to the 
imaginary axis, and from the exterior of the unit circle 
to the positive-real half-plane. The function 



/aH = 



0^ + A2 + 7^(w + xi)'^ + A2 



1 + 7^(1 -a;i)2 + A2 + 7^(1 + m)2 + A2 

(2) 

adds two bumps at symmetric points x'i and —x'i. The 
denominator in Eq. ||2Jl forces /a(w) to map ui = 1 to 1, 
so that fA.e{w) maps infinity to infinity. The parameter 
A controls the size of the bump, and 7 controls the aspect 
ratio of the bump. Since g{w) maps —1 to infinity, we 
choose 6 at random between and tt; and ioi 6 > ^, 
fA.eiw) becomes 



(3) 



where the bar denotes complex conjugation. 

The actual bump positions are not at ±xi, but are off 
by a small factor determined by x and A. The bump size 
is also dependent on x and A. In order to get bumps at 
angle 9, /a must place bumps at g{e^^^) = jf^^i- We 
do this by an approximation method. To keep all the 
particles in the cluster the same size, the bump size on 
the unit circle is varied according to the first derivative 
of the composite conformal map in the original version 
of the conformal map technique fl^ This assumes that 



the higher order derivatives are negHgible, which is not 
true deep inside a fjord. Thus particles added in a fjord 
can end up being very large, and sometimes can partially 
fill the channel. To combat this effect, we measure the 
bump area at each step, and iteratively correct the size 
parameter A if the area is outside of a preset tolerance 
(10% for the results in this paper); compare [l^. 

While conformal mapping allows one to grow DBM 
for any 77, and directly produces a conformal map for 
the cluster boundary, it is computationally intensive. 
A more efficient numerical algorithm for generating off- 
lattice DLA (that is, = 1) is a simple adaptation of 
hierarchical maps |l5l | to channel geometry. This method 
enables close to linear dependence of computing resources 
on cluster size. We used in total 1.7x 10^^ particles for the 
dimension calculations, and (including probes) 4 x 10^^ 
particles for the average profile. 

Both periodic and reflective boundary conditions have 
been implemented on the sides of the channel (the pe- 
riodic boundary condition is sometimes referred to as 
"cylindrical"). The reflective boundary conditions is 
achieved as above: the cluster is grown in a channel of 
double width and periodic boundary conditions, and for 
each deposited particle we deposit also its mirror image. 
At the end one of the images was discarded. By the con- 
ventions used in this paper the channel is given by the 
range — w/2 < y < w/2 and the clusters grow (macro- 
scopically) in the positive x direction. 

Fractal dimension of channel DLA. The fractal dimen- 
sion is measured through the density. The average den- 
sity scales with the width w of the channel, with exponent 
given by the co-dimension: 



p{w) 



,.D-2 



(4) 



To avoid transients, we discarded the first and last 3w 
long section of the clusters, and measured the density 
(number of particle centers per area) on the remaining 
middle section. 

We generated clusters of 8 x 10® to 32 x 10® particles in 
channels of width w = 50, 100, 200, 500, 1000, 2000, 5000 
particle diameters. For each width the number of clusters 
grown ranged from a few hundred to a few thousand, with 
more and larger clusters necessary for large widths, to 
achieve comparable statistical confidence in the average 
density. 

Figure ^ shows the width (w) dependence of the ef- 
fective fractal dimension Z?off = 2 + din p/dlnw. The 
fractal dimension tends to D — 1.712 ± 0.002, indepen- 
dent of the choice of boundary conditions. Off-lattice 
noise reduction does not change the dimension, but 
accelerates the convergence to its asymptotic value. 

Average profile. Analytical solutions for unbranched 
Laplacian growth in a channel with reflective boundary 
conditions have been known for a long time. One can 
find solutions which translate a fixed profile along the 
channel in time. In the absence of surface tension, these 
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FIG. 1: a) The (effective) fractal dimension as a function 
of channel width w. Circles correspond to the original model 
with periodic or reflecting boundary conditions, the triangles 
are made with off- lattice noise reduction [l^. For the three 
cases the ensemble size (and therefore the statistical uncer- 
tainty) was different. The curves correspond to the fitted lines 
in the b) panel, b) The finite size scaling plot of the same 
quantities, see Ref. 17]. Inset: All data is consistent in the 
1(7 — > 00 extrapolation with the dimension D = 1.712 ±0.002. 
(In the inset the data is shifted horizontally for clarity.) 



solutions form a one-parameter family, called Saffman- 
Taylor (ST) fingers 10] . They are parametrized by the 
asymptotic ratio A of the widths of the finger (wfingcr) 
and the channel (w), and have the profile: 



w(l-A) 



2n 



log 



(5) 



Of these solutions A = 1/2 is the most important, because 
in related experiments [li| this profile has been observed 
in the limit of vanishing surface tension. Analytical cal- 
culations jl^i, llSk i20i] show that surface tension — a singu- 
lar perturbation — selects a discrete set of finger solutions 
(only one of which is linearly stable), which all converge 
to the A = 1/2 ST-finger in the limit of zero surface ten- 
sion. 

It has been suggested _21| that the A = 1/2 ST-finger 
solution also models the average profile both of the un- 
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stable (highly branched) Hele-Shaw fingering and of the 
DLA growth in a channel with the corresponding reflec- 
tive boundary conditions. The profile was defined as a 
level set of the ensemble averaged mass density, and for 
the experimental Hele-Shaw profiles half the maximum 
level was used. For DLA growth it was later shown 
that the level set at 0.5 maximum matches the width of 
A = 0.56, whilst the best match to the A = 1/2 profile 
came from the level set at 0.6 maximum. Outside that 
range the authors of Ref. [i^l concluded they could match 
level sets only to finger widths but not to the full shape 
of any ST-finger. 

Here we use a different kind of finger averaging, which 
does not have any fitting parameter (eg. height of level 
set), as follows. In the conformal map method, we di- 
rectly average the map. That is, we choose a set of points 
on the unit circle in the mathematical plane, and repeat- 
edly map to the physical plane. The position of the image 
points averaged over the different maps, that is, over the 
different clusters that we have generated, is a reasonable 
alternative to the ensemble average of |2l|, |22| . 

In Fig. Eta) we show the average conformal map gen- 
erated this way. We see that the average map for DLA, 
77 = 1, does not correspond to the ST result for A — 1/2, 
but we get a good match to it for the DBM growth 
at 77 = 1.2. This is an interesting result, especially in 
the context of recently proposed equivalences between 
DBM models with generalized local spatial cutoff. In that 
framework ^ a higly ramified viscous finger with sim- 
ple surface tension cutoff corresponds to standard (fixed 
size cutoff) DBM with 77 w 1.2. Here we observe that the 
non-branching ST-finger solution is very similar to the 
conformal average of DBM clusters of the same 77. 

To use DLA grown with random walking particles in a 
channel, we need only construct the conformal map from 
the complex unit circle to the perimeter of each cluster, 
and take the average of these maps, as above. The confor- 
mal map is obtained numerically by the following method 
^3- We send M probe particles to the frozen cluster, 
record their impact position, and discard them. These 
points correspond to M uniformly distributed points on 
the unit circle. The landing positions of the probe par- 
ticles are labeled topologically as one encounters them 
when tracking the perimeter of the cluster. Finally the 
m-th point is assigned to the angle 27rm/M of the unit 
circle. This angle has an error of the order of A/^^/^, 
which vanishes for large M. 

We measured the average conformal map on channels 
with reflective boundary conditions, widths ranging from 
10 to 2000 particle diameters. For each width, we grew 
10^ short clusters (only about lOw long), and probed each 
with 10^ test particles. In addition, for a few selected 
widths we probed 10^ clusters with 10® probes each. 

The average map for a wide channel generated this 
way is shown in Fig. [2tb) . The curvature of the tip is 
consistent with that of the A = 1/2 ST-finger (the mea- 
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FIG. 2: a) The average conformal map, generated with it- 
erated conformal maps for 77 — 1, 1.2 and 1.5. The profile 
for rj = 1.2 comes the closest to the ST-finger solution for 
A = 1/2. b) The average map of DLA clusters grown with 
random walking particles in a 1000 particle-diameter wide re- 
flective channel. The profile does not follow any ST-finger 
solution. 



sured curvature of the DLA profile for w = 1000 or 2000 
corresponds to A = 0.51 ± 0.03). The asymptotic width, 
however, is larger. The average map significantly differs 
also from the ST-finger of matching asymptotic width. 

The average conformal map, rescaled onto a unit wide 
channel, shows strong dependence on the channel width. 
This is shown on Fig.|31as a function of reduced x position 
relative to the tip, S. — {x ~ Xtip)/w. Details of the tip 
and tail regions show that these are clearly not consistent 
with a common asymptotic ST-finger shape. 

We performed a finite size scaling on the w depen- 
dence of the finger width. The filling ratio of the finger 

^finger 

/w = [y+{x) - y-{x)]/w was measured at selected 
^ values: ^ = —0.5, —1, and —1.5, and is plotted on 
Fig. 0] The finite size scaling exponent was found to be 
1/ — 0.33, same as for other quantities The most 

interesting is the second extrapolation: ^ — > — cx). We 
only have 3 points for this, so it is only reasonable to 
give a lower bound: u'finger 

/w > 0.62. 

In summary, using large scale random walker based 
simulations we have shown that the fractal dimension 
of DLA in a channel — with either periodic or reflective 
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FIG. 3: The w dependence of the average conformal map, 
rescaled onto a unit wide channel. The tip and the tail are 
magnified on the bottom panels. 
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FIG. 4: Finite size scaling of the finger filling ratio 
'i'flngcr(C)/w'i measured at a few fixed ^ values. Only channels 
50 < w < 2000 were included in the fits. The finite size scaling 
gives the extrapolation w — + oo. For the second extrapolation, 
^ —CO, we can only state that wanga/w > 0.62. 



boundary conditions — is the same as in radial geometry. 
This is a great simplification compared to earUer claims 
of boundary condition (geometry) dependent fractal di- 
mension. Second, using both iterated conformal maps 
and random walker based simulations, we measured the 
average profile of the clusters, defined by the average con- 
formal map, and compared them to ST-finger solutions 
of the corresponding continuum problem. The averaged 
DLA profile is reminiscent but distinct from the ST- 
fingers, while the average profile of DBM clusters with 
rj = 1.2 are rather similar to the ST-finger with A = 1/2. 

We are indebted to Dave Kessler for the suggestion of 
measuring the average conformal map. This research has 
been supported by the EC under Contract No. HPMF- 
CT-2000-00800. The computing facilities were provided 
by the Centre for Scientific Computing of the University 
of Warwick, with support from the JREI. 
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